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An explicit formula for tropospheric range error is derived for a troposphere with a 
spherically symmetric refractivity which varies exponentially with height. The correction 
is given as function of surface refractivity and ray elevation angle, and its accuracy and 
limitations are discussed. 



In the long-range tracking of an object in space by 
an observer on the earth, systematic errors are 
introduced by the earth's atmosphere, most of the 
contribution arising from the troposphere at micro- 
wave frequencies. This paper derives an explicit 
analytic expression for the tropospheric range error 
for a standard atmosphere whose refractivity is 
exponentially distributed with height. 

If one measures the delay 



-H. 



nds 



incurred by a wave between point (0), on the earth's 
surface, and point (1), above the troposphere, then 
the measured range, R e , between (0) and (1) is de- 
fined by 



= nds, 



R e =CT 



the integral extending over the curved ray path. 
Figure 1 illustrates the geometry involved. Here c 
is the velocity of light in vacuo, and n is the index of 
refraction. The true or slant range between (0) and 
(1) is 



#o= J ds, 



the integral being taken along the geometric line 
joining (0) and (1). The error arising from the 
refractivity is 



AR= f 1 nds-Ro. 



(1) 



The range error, AR, arises from (a) the increase in 
path length due to the curvature of the ray path, 
and (b) the decreased propagation velocity. 

It can be shown that the error arising from the 
increase in path length is of the second order in 
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('ft— 1), whereas the error arising from the decrease in 
wave velocity is of the first order in (n-l), so that as 
a first approximation, 
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A 1 R= C (n-l)ds, 


(2) 


where the integral is now taken along the slant range 
path. This mathematical stratagem permits one to 
integrate the value of AJI directly when n has a 
simple functional dependence on height. The valid- 
ity of this approximation is evident from table 1 
which enumerates the range error, AR (column 2), 
the excess of ray path over geometric path, Ai 
(column 3) and percent range error due to the excess 
of path length, 100 AL/AR (column 4), as a function 
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Figure 1. Geometry of ray p 


ath. 



of ray initial elevation angle, O . These values were 
computed from ray tracings for a typical surface 
refractivity N s =(n s — 1)X10 6 = 313 for an exponen- 
tial reference atmosphere and are taken from an 
article [1959a] by Bean and Thayer. 

Table 1 shows that even for the worst possible case, 
only 8 percent of the total range error arises from the 
excess of the ray over the geometric path length, and 
that this percentage decreases rapidly as the initial 
ray elevation angle increases, becoming negligible for 
initial ray elevations of a few degrees. 

Table 1. 



0o, Initial 


A .R= range 


Excess of ray 




ray 


error 


path over 


100 A LIAR 


elevation 




geometric path 




mr 


m 


m 







104 


9 


8.2 


8 


81.4 


4.8 


5.8 


15 


68.1 


3.2 


4.4 


30 


49.7 


1.3 


2.6 


65 


29.5 


0.3 


1.0 


100 


20.7 


0.1 


0.5 


200 


10.9 






400 


5.6 







AR islhot very sensitive to the distribution of 
n-l=NX10~ 6 with height, since only the integrated 
effect or the first moment of N comes into considera- 
tion. Also, it has been found [Bean and Thayer, 
1959b; Hoehndorf, 1961] that the average decay of 
N with height is nearly exponential, so that the dis- 
tribution of N may be approximately described by 



N=N s exp {-c(h-h s )}, 



(3) 



where N s is the surface refractivity at a locality 
which is height h s above sea level. Further exten- 
sive analyses [Bean and Thayer, 1951 a and b] have 
shown that there is a very strong correlation between 
N s , AN, and h s of (3), where AiVis the decrease from 
surface refractivity observed 1 km above the surface. 
Utilizing this mean relationship between AiV andiVs, 
Bean and Thayer have published tables giving c as a 
function of N s . Other methods of determining the 
value of the parameter c which will best enable one to 
predict the tropospheric range error have been ex- 
plored. However, if one assumes that the distribu- 
tion of N is exponential with height, and if one uti- 
lizes the approximation given in (2), then (2) may 
be integrated explicitly as is shown in the appendix, 
giving 
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Table 2. Comparison o/AiR and AR for N s = 313 for target 
above troposphere 
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and 



#=tan p \cr /2. 



(5) 



Here r is the distance from earth's center to ground, 
that is, the earth's radius plus height above sea level, 
£0 is the slant range elevation angle at ground. 

For the point (1) situated above the troposphere 
the range of integration of (4) may be extended to 
infinity with negligible error, yielding 



A^IO" 6 Nds= — r y [1-erf (#)]. 
Jo c sin p 



(6) 



For g large, (4) approaches the plane earth approx- 
imation, 



A 2 R= 



lO-Wo 



c sin ft 



(7) 



which, of course, does not converge at ift=0, whereas 
(6) does. 

Table 2 enumerates the different computed range 
corrections, where Aii? is the approximate correction 
(6), Ait the true correction, and the last column is 
the plane earth correction. 

Comparing A L 7? and AR, it is apparent that (6) is 
a good approximation to the tropospheric error one 
would compute for an exponential atmosphere from 
a ray tracing. Its simple form suggests that it will 
be valuable for predicting tropospheric range error 
from a knowledge of only the surface conditions. 
The question of how to determine the decay constant, 
c, to give accurate predictions and the probable error 
of such predictions will be reported in a future paper. 

It is worth remarking that for 

N=N s exp {-c(h-h s )-d(h-h s ) 2 } 

the expression for range error, A^, turns out to be 
precisely (4) and (6), except that now one has 

c sin ft 

V2c cos 2 0o , A i • 2 o 
+4dsin 2 ft 
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Appendix 

Consider finds. Since the difference between the 
distance along the ray path and that along the slant 
line path is negligible (see fig. 1), 

r<nds=r- f^s+lO" 6 rNds=R +10- Q ^Nds, (Al) 
Jo Jo Jo Jo 

where n=l+iVX10~ 6 , R is the slant range, and the 
integrals on the right hand side are taken along the 
slant range path. Along this geometric line, 



d<p 



r cos O 



dr 



r-\Jr 2 —rl cos 2 f$ 



(A2) 



where r, p are polar coordinates with earth's center 
as origin, and designates the angle between slant 
path direction and the tangent to the radial, that is, 
the elevation angle. r , O pertain to the initial path 
point. Then, letting K=r cosl3 , 



8 (r) = J^^-J*=K* 



(A3) 



is the distance along the geometric path. Expressing 
r in terms of s, and letting h'=h—h s designate 
height above the earth's surface, 

ii/ ["-, > s 2 , 2s sin0 o "]2 
r=r +h'=r Q [l+^+ ^ ° J ■ (A4) 



Expanding in terms of s/r to terms of the second 
order, 



W = s sin P +^- cos 2 o . 
zr 



(A5) 



Letting N=N s e- ch ', 
f 1 Nds=N s r° exp / -cfs 2 ^-^+s sin O ~|\ ds. 



(A6) 



Integrating (A6), 



-t*dt 



(A7) 



(A8) 



where erf (#)h=— e~ { ' 

■\'7T Jo 

and where 9—^J^f tan £ . 
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